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Abstract. We show how easy it is to get the Sommerfeld enhancement for a Yukawa potential, 
for definite partial waves, beyond the S wave analyzed in previous literature. In particular, we report 
results for the P wave (for which there is a resonant pattern and the enhancement can be of several 
orders of magnitude even far from resonances) that could be relevant for the analysis of experimental 
cosmic rays data possibly signaling the annihilation of dark matter particles. 

^^ ■ 

^ ! 1 The general formula 

\ Here we derive the formula for the Sommerfeld enhancement for a definite partial wave / in the case 

■ of an attractive Yukawa potential — -^^^r— • 

' This is the non-relativistic limit of a process in which two particles attract each other by ex- 

changing repeatedly a massive vector (or also scalar) boson, before undergoing some annihilation 
reaction. The Sommerfeld enhancement could be relevant for the interpretation of cosmic rays 
data as due to the annihilation of dark matter particles, and it has been invoked in various pa- 
pers [1], [2], [4], [3], [5], [6], [7], [8]. However, as far as we know, it has not been studied for partial 
^ , waves higher than I = 0, in particular for the I = 1 V partial wave, which could be particularly 

relevant. 

In general one gets a matrix valued potential. In the analysis below we assume that that the 
pH ' matrix can be diagonalized and that it has a Yukawa form for the eigenstates of the two particles. 

I— We consider therefore a definite eigenstate. 
' The starting point is eq.(2.35) of ref [9]: 

m ■ 

<^ ' where Ai{p,p') is the I partial wave amphtude for the process described above, p being the CM 

^f-^ . momentum of the interacting particles and p' the one of the reaction's products; 

' Aq^i{p,p') = p^aQj{p') is the bare (i.e. for a = 0) partial wave reaction amplitude. 

0^ . We have assumed that the bare reaction amplitude is proportional to p', which, for small p, is 

' the typical dependence on the incoming momentum of the / partial wave. 

^ , Rp^i{r) is the solution of the partial wave Schrocdinger equation: 

1 fRp^ 2dRp,i l{l + l)R^ p^ , „ 
- 2^^^^;;^ r—r 7^ ^ " W ^ " ^ ^ 

normalized such that 

poo 

r^drRg,i{r)Rp,iir)^d{p-q) (3) 
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Here = rn/2 is the reduced mass and m is the mass of the two particles. 

We know from textbooks (see for instance [10]) that this normalization corresponds to the asymp- 
totic behavior _ 

2 sm{pr - if + Si) 



RpAr\^oo^\l-—^ (4) 
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Let us define x ~ pr and put Rpj{r) ~ Npx^i{x); the equation for <I>; is: 

^ ^([±11 + + 1)$^ = (5) 

where a = a/v, b = fi/ {rrirv) and v = p/rur is the relative velocity. 
Suppose we solve this equation with the initial conditions 

$K0)-1 <i>;(0) = -a/(/ + l) (6) 

(the condition for <1>J(0) is dictated by the equation for a regular solution). Then the asymptotic 
behavior will be 

x'+''^i{x),^^^CsHx-^^5i) (7) 

In order to agree with the normalization of eq. Q we have to put N = \J^-^ ■ Substituting in eq. ([T]) 
wc get 

Mp.p) = ^ ^p'ao,i = ^ '-Ao,i{p,p') (8) 

In conclusion, by defining the Sommcrfeld enhancement enhi for the I partial wave cross-section (or 
equivalently for the rate) as 

ai = enhi ■ ao,i (9) 

wc get 

, l-3---(2/ + l) ^, 
enhi = ( ^ ) (10) 

where C is obtained by looking at the asymptotic behavior eq.© of the solution of eq.® with 
the initial conditions eq.®. enhi depends on the two parameters a and b. It is not necessary to 
determine 6i. 

Another equivalent strategy is to put Rpj{r) — Npipi{x)/x; the equation for ipi is: 

^r + (l + -e-^-^V.-0 (11) 
X x^ 

If one solves this equation with the initial conditions corresponding to 

fi{x)^^,~.x'+^ (12) 

then the asymptotic behavior will be 

111 

^i{x).^^^^Cfi\ii{x- — +5i) (13) 
with the same C of eq.® giving the enhancement as in eq. pop . 



2 Computations for / = 1 

In principle it is easy to get C: for instance one can use the NDSolve instruction of Mathematica to 
get the numerical solution of eq.® with the initial conditions eq.®, or equivalently of eq. pip with 
initial conditions eq. . In order to find C one takes Fi{x) = a;'+^$/(a;) or else F{x) = (p{x), and 
one plots Fi{xY + Fi{x — 7r/2)^ for large x: when this is constant it is equal to . We follow the 
strategy of eas. (|llll2p . which provides more clean numerical results. 
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Figure 1: Fo{x)'^ + Fa{x - ixj^f for a = 10"^ a* = IGei', m = lO^Gew, v = 2 x 10^ 
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Figure 2: F\{x)^ + — (left) and ;/(.t) (right) for the same parameters of Fig.l. 

It is expected that this procedure should work less well for h very low and a very large because in 
this case the asymptotia is reached for very large x and the numerical solution accumulates errors. 
However, for 6 = and for any a we already have the exact result derived analytically in ref. [9]: 



sinh{'Ka)l 
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(14) 



In practice, this works well for I = 0, see for instance Fig. 1. 

For I = 1 the quantity F/(a;)^ + Fi{x ~ 7r/2)^ may sometimes continue to show decreasing oscil- 
lations: here it is convenient to take into account the sub-leading term in the asymptotic expansion 
(fT5|l . which we know to be the free wave function up to the phase shift 6i. Therefore, the improved 
version of (fT3l) for / = 1 is: 



C ■ ( sin(a; 



cos X 



{x-'^ + 5i). 



(15) 



To get C^, one defines, in three steps, with Fi{x) :~ tpi[x), 



k{x) 
h{x) 



(tt^ - 16x2) 



2\2 



-{Fi{x + Ti/Af + Fi{x - 7r/4)2), 



(87r(7r2 - 16x2))' 
(tt^ - 4:TTX^){k{x + 7r/4) + k{x - 7r/4), 
h{x + n/A) + h{x ~ Tr/4) 
(7r2-16x2)(37r2 + 16(l + x2))' 



For X large (say x > 30), j{x) quickly converges to a constant which equals C^. (For I > 1 the free 
wave function is more complicated and one should do more steps). It must be said that the results 
of the improved procedure differ little from what could be obtained simply by finding by eye the 
average of the oscillations of Fi{x)'^ + Fi{x — 7r/2)2, see Fig. 2. 
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Figure 4: enhi as a function of m{Gev) for Vsingie particle = 10 ^ blue, 10 ^ red , 10 ^ green. Here 
a = 1/100, n = IGev. 



The case of the S / = wave has been discussed at length in the literature, and it has been found 
a resonant pattern, see in particular refs. [6] and [8]. We have verified that we get the same results. 

Here we present some numerical result for the P / = 1 wave, which also shows a resonant pattern. 

In Fig. 3 we show the enhancement ea. pH)) for / = 1, taking the values of the parameters used in 
the numerical evaluations for the S wave in ref. [8] for the range of ref. [6], that is a = 1/30, fi = 
90Gev, as a function of m (expressed in Gev) for Vsingie particle = 10"'^, 10~^, 10~^. 

In Fig. 4 we show the enhancement eq. ([T(I|) for ^ = 1, taking the values of the parameters used in 
the numerical evaluations for the S wave in ref. [8] for the range of ref. [5], that is a = 1/100, ^ = 
IGev, as a function of ra (expressed in Gev) for v^ingie particle = 10"'^, 10~*, 10~^. 
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